Four-fold basal plane anisotropy of the nonlocal magnetization of YNi2B2C 
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Studies of single crystal YNi2B2C have revealed a four-fold anisotropy of the equilibrium magne- 
tization in the square crystallographic basal plane. This 7r/2 periodicity occurs deep in the super- 
conductive mixed state. In this crystal symmetry, an ordinary superconductive mass anisotropy (as 
in usual London theory) allows only a constant, isotropic response. In contrast, the experimental 
results are well described by generalized London theory incorporating non-local electrodynamics, as 
needed for this clean, intermediate-^ superconductor. 



Borocarbide superconductors have received consider- 
able recent attention, due in part to the interaction be- 
tween magnetism and superconductivity. A rich super- 
conducting phase diagram, including transitions between 
hexagonal, rhombohedral, and square vortex lattices has 
been observcdEl^Q. The existence of vortex lattices with 
non-hexagonal symmetry has been attributed to nonJa- 
cality effects on the superconducting electrodynamicsQfl, 
which arise from the large electronic mean free path of 
these clean superconductors. Geometrically, a vortex di- 
rected along the tetiagonal c-axis has cuixent contours 
that are square-likaJ. It has been shownQ in the non- 
magnetic borocarbide YNi2B2C, that the deviations from 
the standard (local) London magnetic field dependence of 
the equilibrium magnetization Meg oc ln(i/) can be quan- 
titatively accounted for by introducing aon-local electro- 
dynamics into the London descriptionEl. Traditionally, 
it was widely thought that nonlocality effects should be 
significant only in materials with a Ginzburg-Landau pa- 
rameter K ~ 1; however, YNi2B2G has k « 10 - 15. 

In the local London model of superconducting vor- 
tices, the material anisotropy is introduced via a second 
rank mass tensor m^ . In tetragonal materials such as 
YNi2B2C or LuNi2B2C, the masses in both principal di- 
rections in the square basal plane are the same, rua = rrib, 
thus the superconducting properties are isotropic in the 
a-b plane. In-contrast, non-local corrections are expected 
to introduceO a four-fold anisotropy as a function of the 
magnetic field orientation within the a-b plane. A tem- 
perature dependent in-plane anisetropy of the upper crit- 
ical field Hc2 has been observedll2l in the non-magnetic 
borocarbide LuNi2B2C and described within a Ginzburg- 
Landau framework incorporating non-local effects. How- 
ever, a direct observation of the in-plane anisotropy deep 
into the superconducting phase, where the non-local Lon- 
don model applies and unusual vortex lattices are ob- 
served, has not been reported until now. 

In this work we show that, in the superconducting 
mixed state of YNi2B2C, the reversible magnetization 
oscillates with a 7r/2 periodicity when the applied field is 



rotated within the a-b plane. The amplitude of the an- 
gular oscillation decreases with field, passes through zero 
and then reverses sign at a field well below Hc2- The 
results are in good quantitative agreement with the wn- 
local London description introduced by Kogan et al.H 

The 17 mg single crystal of YNi2B2C investigated in 
this study is the same as that previously used by Song et 
al.Ll to explore the magnetic response when the applied 
field H is parallel to the c-axis. The critical temperature 
is Tc = 14.5-?^, defined as the point at which the linearly 
varying magnetization M{T) extrapolates to zero; this 
ignores a slight "tail" extending to Ih.QK . The crystal is 
a slab of thickness t ^ O.bmm, whose shape and size in 
the basal plane are sketched in Fig. 1. It will be useful 
to approximate such shape by an ellipse of axis and 
Ly. X-ray diffraction shows that the crystallographic c- 
axis is normal to the slab, and that the two equivalent 
(110) axes of the tetragonal structure very approximately 
coincide with the axes of the ellipse. 

Measurements were performed in a Quantum Design 
SQUID magnetometer with a 50 kOe magnet. Two sets 
of detection coils allows us to measure both the longi- 
tudinal (parallel to H) and transverse (perpendicular to 
H) components of the magnetization vector M, but only 
the longitudinal component (denoted hereafter as simply 
the magnetization M) will be discussed in this vprk. The 
crystal was mounted into a previously describedLiJ, home- 
made rotating sample holder with rotation axis perpen- 
dicular to H. The c-axis was aligned with the rotation 
axis, so that H could be rotated within the basal plane. 

Figure 1 shows M as a function of the angle (p between 
H and the a-axis (see sketch in Fig. 1), at T — 7K 
and two values of H. The crystal was initially cooled 
in zero field, H was then applied and the sample was 
subsequently rotated in steps of Aip « 3.1°. Each data 
point was taken at fixed ip. 

We first discuss the low field curve of fig. lb. As H = 
30Oe is well below the lower critical field Hd for all (p, 
this curve represents the total flux exclusion of the Meiss- 
ner state. The oscillatory behavior with periodicity tt 
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(two-fold symmetry) originates from geometrical effects. 
Indeed, a field applied at any orientation within the basal 
plane can be decomposed in = Hcos{if + 45°) and 
Hy — Hsin{ip+4:5°). If we approximate the crystal shape 
by the ellipse, the Meissner response associated with each 
component is At: Mi = —Hi/ {I — Vi), where i = x;y and 
i^i = t/Li are the demagnetizing factors, thus AttM — 
-H[cos'^{(p + 'i5°)/{l-i^^) + sin^{ip + A5°)/{l-iyy)]. The 
best fit to this expression gives I'x ^ 1/4 and i^y ^ 1/5. 
This corresponds to the ellipse of axes ~ 2.0mm and 
Ly 2.5mm shown in the sketch of fig. 1. 

We now turn to the high field data of Fig. la. The 
applied field, H = AbkOe, is well above Hc2 ~ 35 kOe at 
this temperature (which is only weakly (p dependent, see 
below), thus in this case M{ip) = M"*'((^) is the normal 
state paramagnetic response. We again observe an oscil- 
latory behavior, but in this case the periodicity is 7r/2. 
By combining the information provided by the X-rays 
with the geometrical effects on the Meissner response, 
we conclude that the maximum normal state magneti- 
zation occurs at the crystallographic orientations (110) 
and (110), while the minimum corresponds to (100) and 
(010). No hint of the geometry-originated two- fold sym- 
metry is observed. This is to be expected, as demag- 
netizing effects vanish in the limit \M/H\ ^ 1. Further 
analysis of Af"'* suggests that it arises from a low concen- 
tration, ~ 0.001 molar fraction, of rare earth ions, most 
likely from impurities in the yttrium starting metal. As 
shown below, Af"** is much smaller than the supercon- 
ducting contribution except close to Hc2- 

The above procedure cannot be used in the supercon- 
ducting mixed state, due to the appearance of magnetic 
hysteresis arising from vortex pinning. The iritical cur- 
rent density Jc is very small in this crystalEI. As a re- 
sult, the magnetic hysteresis {M^ — AfT) ex; Jc, where 
and Af^ are respectively the magnetizations measured in 
the field-decreasing and field-increasing branches of an 
isothermal M{H) loop, is small as compared to the equi- 
librium or reversible magnetization, Afeg « {M^+M^)/2. 
In spite of this, the residual hysteresis strongly affects the 
response obtained by rotating the crystal at fixed T and 
H, by superimposing a periodicity tt (related to shape 
effects on the critical state magnetization) that almost 
completely hides the intrinsic 7r/2 periodicity of funda- 
mental interest. 

To solve this difficulty, we performed magnetization 
loops at T = IK at a set of fixed angles and then cal- 
culated Meq{H) for each (p. In all cases we extended 
the loops up to iJ = 50fcOe, thus we could repeat the 
measurement of Af"'* in the normal state and compare 
the data with those obtained by rotating at fixed H. 
Due to the absence of hysteresis, both determinations 
of M^^^H, (p) should coincide. This is indeed the case as 
seen in Fig. la, where the open circles represent the data 
a.t H = 45kOe obtained from the M{H) loops. 

Figure 2 shows Mf.q (obtained from averaging M'^ and 
M'^) as a function of ip for several H. All the data have 
the same scale, but the curves at different H have been 



vertically shifted to accommodate the whole field range 
within the plot. For H < l.bkOe, the irreversibility 
becomes large enough to introduce a significant uncer- 
tainty in the determination of M^q, consequently those 
data have been disregarded. It is apparent that a four- 
fold symmetry exists in the whole field range of the 
measurements. To quantify the amplitude of the oscil- 
lations, we fitted the curves by Meq{H,(p) = (Meq) + 
5Mcq{H)cos{Aip). 

The values of 6Meq{H) so obtained are plotted in Fig- 
ure 3, while the values of {Meq){H) are shown in the 
inset. A remarkable fact, clearly visible in figs. 2 and 3, 
is that 6Meq{H) crosses zero and it reverses sign at some 
intermediate field {r^ 12 kOe) well within the supercon- 
ducting mixed state. Another interesting observation is 
that the amplitude of the oscillations at 7J ~ 1.5 - 2 kOe 
is as large as that ai H ^ 50kOe. 

The above results show that a 7r/2 basal plane 
anisotropy exists both in the normal and superconduct- 
ing states. It is also clear from Fig. 3 that a change 
in the behavior of SMeq{H) takes place at the supercon- 
ducting transition at Hc2 ^ 35fcOe. This observation, 
together with the sign reversal and the large amplitudes 
at low fields, point to the existence of a second source of 
in-plane anisotropy, in addition to the normal state one, 
that turns on in the superconducting phase. 

Prior to analyzing the superconducting basal-plane 
anisotropy it is necessary to subtract the normal state 
contribution, which persists within the superconducting 
phase. To that end we performed rotations at fixed as 
those shown in Fig. la, at several T and H above Hc2{T). 
We found a paramagnetic response that exhibits a four- 
fold symmetry, with the minimum at = in all cases, 
i.e., Af"^((/3 = 45°) > A4^"*(^ = 0°) > for all T and 
H. We thus have a well defined set of data SM"'{H,T) 
which exhibits no sign reversal. The extrapolation is not 
obvious, however, as SM^" is not linear in H . Figure 4 
shows all the ^Af"^ data collected at various tempera- 
tures 5K < T < 16K and H < 50fcOe, as a function of 
H/T. We found that, when plotted in this way, all the 
data points collapse on a single curve. 

The dashed line in Fig. 4 is a fit to the SM"'{H/T) 
data. The same fit, for the case T — 7K, is also shown as 
a dashed line in Fig. 3. We can now subtract that curve 
from the total SM^q shown in Fig. 3, to isolate the su- 
perconducting contribution SM^^. Note that, as JAf""* is 
always positive and increases monotonically with H, both 
the sign reversal and the non-monotonic behavior imme- 
diately below Hc2 exhibited by dMeq must arise from the 
SM^q contribution. 

We now show that the four- fold symmetry of M^q, as 
well as the field dependence of SM^^, can be well de- 
scribed using the non-local modifications to. the London 
electrodynamics introduced by Kogan et al.B. According 
to that model, for Hd <C Hc2 
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Here Mq = $o/327r^A^, the new characteristic field 
Hq = $o/47r^p^ is related to the nonlocality radius p, 
and C — Vi ~ ^n{HQ/rj2Hc2 + 1), where rji and 772 are 
constants of arder unity. 

Song et aLQ have shown that the magnetization of this 
same crystal is very well described by Kogan's model, 
when H \\ c-axis. A fingerprint of the nonlocality effects 
is the deviation from the oc In(iJ) behavior predicted 
by the local London model. The curvature clearly visible 
in the inset of Fig. 3 is thus a strong indication that 
nonlocality also plays a major role when H _L c-axis. It is 
worth mentioning that, although a quantitative analysis 
of the curve in the inset of Fig. 3 in terms of Eq. |l] would 
require removal of the normal state magnetization, its 
contribution is small and would not significantly modify 
the curvature seen in the M^q vs. In(iJ) data. 

We now apply Eq. ^ to the analysis of our data. In 
principle, the basal plane anisotropy could be ascribed 
to the material parameters Mq, i?o and Hc2- However, 
Mq (X is isotropic within the a-h plane of a tetragonal 
structure. On the other hand, four-fold variations of Hc2 
within the basal plane have been observed in LuNi2B2C 
and attributed to nonlocalityll3. We then assume that 
both Hq and Hc2 have 7r/2 periodicity, Ho{(p) = {Hq) + 
SHocos{4:(p) and Hc2{<f) = {Hc2) + 6Hc2Cos{4:ip). To first 
order in SHq and SHc2 we obtain 



fit deviates from the data close to Hc2, where the Lon- 
don model fails. These experimental results show that 
nonlocality effects have a profound effect on these clean, 
intermediate-K superconductors and they underscore the 
remarkable utility of the generalized London theory. 

In summary, we have demonstrated a four-fold 
anisotropy in the square basal plane of clean single crys- 
tal YNi2B2C. This superconducting response is incon- 
sistent with conventional local London theory, but it is 
well explained by a generalized London model incorpo- 
rating non-local electrodynamics, with parameters based 
largely on complementary experiments. These observed 
effects of nonlocality persist deep into the superconduct- 
ing state, where complex, evolving vortex lattices occur. 

[Note added: while preparing this manuscript, we 
learned that P. C. Canfield et al. at the Ames Lab. have 
observed similar oscillations of the basal plane magneti- 
zation of LuNi2B2C.] 
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Experimentally we have determined {Hc2) — 35 kOe 
and 5Hc2 ^ 0.4 kOe (at T = 7K), so we can fix £2 = 0.01. 
We could also attempt to determine Mq and (Hq) by fit- 
ting our Meq datEuwith Eq. |l|. However, this is a difficult 
task that requiresQ the measurement of a large set of tem- 
peratures to check the consistency of the results. Instead, 
we decided to use the results of Song et aLQ, (for H \\ c- 
axis) as good estimates. For T = 7 K, we take (Hq) = 56 
kOe and Mq = 5.2 G. (Here we scaled down Mq oc 1/XaXc 
by the experimental mass anisotropy, 7 w 1.15, between 
the c-axis and the a6-plane, which is close to tha-value 
7 ~ 1.1 obtained from band structure calculationsta. ) In 
any case, small variations in any of these parameters will 
not significantly affect the rest of the analysis. If we also 
assume 772 1, we obtain a ~ 2/3. With these fixed 
parameters, we fit our SM^q{H) data with Eq. ^, with 
the single free parameter ei. We obtain ci — 0.14. The 
fitted curve (for (5Af|^ -I- SM"'^) is shown as a solid line 
in Fig. 3. 

The very good coincidence between our data and the 
model is remarkable. With a single fitting parameter 
ei, which is field independent, we have been able to ac- 
count for the nontrivial H dependence of SM^^, including 
the sign reversal at intermediate fields. Of course, the 
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FIG. 1. Equilibrium magnetization M^q at T = 7K, as a 
function of the angle between the applied field H (contained 
in the a-b piano) and the a-axis, for (a) H — ASkOe and (b) 
H = 30Oe. The inset shows a sketch of the crystal shape and 
crystallographic orientation in the basal plane, superimposed 
to an ellipse with the same demagnetizing factors. 

FIG. 2. Reversible magnetization M^q at T — 7K, as a 
function of the angle (fi between the applied field H (contained 
in the a-b plane) and the a-axis. The fields (in kOe) are 
indicated next to each curve. The scale is the same for all the 
curves, but data at different H have been vertically shifted. 

FIG. 3. Amplitude SMeq of the four-fold oscillations of the 
basal plane magnetization, as a function of the applied field 

H. The dashed line is the normal state contribution 
obtained from the fit shown in Fig. 4. The solid line is the 
fit to SM;;^ + SM"^" using Eq. (2) with ei as the only fitting 
parameter. Inset: average in-plane magnetization {Meq) as a 
function of applied field (see text). 

FIG. 4. Amplitude SM"" of the four-fold oscillations of 

the basal plane magnetization in the normal state, as a func- 
tion of H/T. The dashed line is a polynomial fit. 
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